Non linear transport theory for negative-differential resistance states of two 
dimensional electron systems in strong magnetic fields. 
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We present a model to describe the nonlinear response to a direct dc current applied to a two- 
dimensional electron system in a strong magnetic field. The model is based on the solution of 
the von Neumann equation incorporating the exact dynamics of two-dimensional damped electrons 
in the presence of arbitrarily strong magnetic and dc electric fields, while the effects of randomly 
distributed impurities are perturbatively added. From the analysis of the differential resistivity and 
the longitudinal voltage we observe the formation of negative differential resistivity states (NDRS) 
that are the precursors of the zero differential resistivity states (ZDRS). The theoretical predictions 
correctly reproduce the main experimental features provided that the inelastic scattering rate obey 
a T 2 temperature dependence, consistent with electron-electron interaction effects. 

PACS numbers: 73.43.Qt,71.70.Di,73.43.Cd,73.50.Bk,73.50.Fq 



I. INTRODUCTION 

In the past few years the study of non-equilibrium 
magneto-transport in high mobility two-dimensional 
electron systems (2DES) has received much attention 
due to the experimental finding of intense oscillations of 
the magneto- resistivity and zero resistance states (ZRS). 
Microwave-induced resistance oscillations (MIRO) were 
discovered 1 -' 2 -' 3 -' 4 - in 2DES samples subjected to microwave 
irradiation and moderate magnetic fields. For the MIRO 
the photoresistance is a function of the ratio e ac = oj/u> c 
where u) and lo c are microwave and cyclotron frequencies. 
This outstanding discovery triggered a great amount of 
theoretical work 5 i 6 i 7 i 8 i 9 i 10 i u i 12 i 13 i 14 i 15 i 16 i 17 . Our current 
understanding of this phenomenon rests upon models 
that predict the existence of negative-resistance states 
(NRS) yielding an instability that rapidly drive the sys- 
tem into a ZRS 1 -^ Two distinct mechanisms for the gen- 
eration of NRS are known, one is based in the microwave- 
induced impurity scattering 5 -' 7 -' 9 - ' 10 ' 11 ' 12 ' 13 , while the sec- 
ond is linked to inelastic processes leading to a non-trivial 
distribution functiona l 14 ' 15 ! 17 . 

An analogous effect, Hall field-induced resistance oscil- 
lations (HIRO) has been observed in high mobility sam- 
ples in response to a dc-current excitatio n 19 ' 20 ' 21 . Al- 
though MIRO and HIRO are basically different phenom- 
ena both rely on the commensurability of the cyclotron 
frequency with a characteristic parameter; in both cases 
oscillations are periodic in l/B. In HIRO the oscillation 
peaks, observed in differential resistance, appear at inte- 
ger values of the dimensionless parameter e dc = ujh/oj. 
Here, Hujh ~ &Eh(2Rc) is the energy associated with 
the Hall voltage drop across the cyclotron diameter; Eh 
is the Hall field and Rc the cyclotron radius of the elec- 
tron at the Fermi level. It has been found that there are 
two main contributions to the HIRO: the inelastic one is 



related to the formation of a non-equilibrium distribution 
function component that oscillates as a function of the 
energy^ 2 - and the elastic contribution is related to elec- 
tron transitions between different LLs due to impurity 
scattering 2 ^. The first one was shown to be dominant at 
relatively weak electric fields, and the latter prevails in 
the strong-field regime. 

More recently it has been demonstrated that the effects 
of a direct dc current on electron transport can be quite 
dramatic leading to zero differential resistance states 
(ZDRS) 2 ^ 5 -. As compared with the HIRO conditions, 
the ZDRS are observed under dc bias at higher magnetic 
fields (0.5- LOT) and lower mobilities (70 - 85 m 2 /Vs). 
At low temperature and above a threshold bias current 
the differential resistivity vanishes and the longitudinal 
dc voltage becomes constant. Positive values for the dif- 
ferential resistance are recovered at higher bias as the 
longitudinal dc voltage slope becomes positive. Bykov 
et al. analyzed the results following an approach similar 
to that of Andreev et ali 1 ^; the presence of the ZDRS is 
attributed to the formation of negative differential resis- 
tance states (NDRS) that yields an instability that drives 
the system into a ZDRS. Similar results where obtained 
by Chen et ali2& 

In this paper we present a model to explain the for- 
mation of NDRS. According to our formalism both the 
effects of elastic impurity scattering as well as those re- 
lated to inelastic processes play an important role. The 
model is based on the solution of the von Neumann equa- 
tion for 2D damped electrons, subjected to arbitrarily 
strong magnetic and dc electric fields, in addition to the 
weak effects of randomly distributed impurities. This 
procedures yields a Kubo formula that includes the non- 
linear response with respect to the dc electric field. Con- 
sidering a current controlled scheme, we obtain a set of 
nonlinear self-consistent relations that allow us to deter- 
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FIG. 1: Differential resistance r xx as a function of the dc bias 
J x for B = 0.784T temperatures from T = IK to T = WK. 
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FIG. 2: Electric field E x as a function of the dc bias J x for 
B = 0.784T and for fixed temperatures ranging from T = IK 
to T = WK. 



mine the longitudinal and Hall electric fields in terms of 
the imposed external current. It is shown that in or- 
der to correctly reproduce the main experimental results 
the inelastic scattering rate must obey a T 2 temperature 
dependence, consistent with electron-electron Coulomb 
interaction as the dominant inelastic process. 



II. MODEL 

We start with the Hamiltonian for an electron in the ef- 
fective mass approximation in two dimensions subject to 
a uniform perpendicular magnetic field B = (0, 0, B), an 



in-plane electric field E = (E x ,E y ,0), and the impurity 
scattering potential V. Hence the dynamics is governed 
by the total Hamiltonian H = H e + V, with 

H e = H Q + eE ■ x , (1) 

here H = II 2 /2m, m is the effective mass of the electron, 
e is the electron's charge, II = p + eA is the velocity 
operator and the vector potential in the symmetric gauge 
is given as A = (—By,Bx) /2. The impurity scattering 
potential is expressed in terms of its Fourier components 

V (r) = e-"l*l / 7-^2 V (?) cxp [iq ■ (r - r,)] , (2) 
i J ( 27r ) 

where r, is the position of the ith impurity and Ni is 
the number of impurities. The explicit form of V (q) de- 
pends on the nature of the scatterersi 2 -, for simplicity we 
assume short-range uncorrelated scatterers. The factor 
exp(— ?]\t\) takes care of the adiabatic switching of the 
impurity potential at the initial time to — > — oo. 

The motion of a planar electron in magnetic and 
electric fields can be decomposed into the guiding cen- 
ter coordinates Q and the relative coordinates R = 
(— H y ,Il x ) /eB, such that the position of the electron is 
given by r = Q + R. The guiding center coordinates is 
written as Q = (Q x , Qy) /eB, in the symmetric gauge 
(Qxi Qy) — (Px + eBy/2,p y — eBx/2). The commuta- 
tion relations for velocity and guiding center operators 
are [n x ,n y ] = [Qn, <2 y ] = —iheB, with all the other 
commutators being zero. 

Our aim now is to compute the electric current density. 
In order to calculate the expectation value of the current 
density we need the time-dependent matrix p(t) which 
obeys the von Neumman's equation ihdp/dt = [H,p\. 
We assume that in the absence of the impurity potential 
the density matrix reduces to the equilibrium density ma- 
trix given by po = f(Ho), with f(E) given by the Fermi 
distribution function. In order to solve the von Neum- 
man's equation we apply three unitary transformations: 
the first two transformations exactly take into account 
the effects of the electric and magnetic fields, whereas 
the third transformation incorporates the impurity scat- 
tering effects to second order in time dependent perturba- 
tion theory. First we consider the unitary transformation 

W(t)=eif Cdt e- l ^-e l -^e t ne l ^ 3L (3) 

where v x (t) , v y (t) , X (t) and Y (t) are solutions of the 
dynamical equations 

1 e -Eii 

v x + -v x + u c v y + —E x =0, X--£ = 0, (4) 
Tj m B 

1 e -Ex 

i' y + —Vy - uj c v x + — E y = 0, Y + — = 0. (5) 

Except for the damping terms, these equations follow 
from the variation of the classical Lagrangian C 12 . The 
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variables v x and v y correspond to the electron veloc- 
ity components and X and Y are the coordinates that 
follow the drift of the electron's orbit. In order to in- 
corporate dissipative effect we added the damping term 
v/ri the dynamical equations. This procedure yields a 
simple scheme to incorporate dissipation to the quan- 
tum system. Recent magnetoresistance experiment a 27 i 28 
and theory 2 ^ suggest, that in 2DES, electron-electron in- 
teraction provide an important contribution to the in- 
elastic scattering rate, giving rise to 1/r, oc T 2 tem- 
perature dependance. Consequently, in what follows we 
shall assume that the inelastic scattering rate is given 
by l/n re (fc B T) 2 /^j4aaa, where E F is the Fermi 
energy. 

The transformation (|3|) renders von Neumann equation 
into the following form 

The electric field term is conveniently removed from the 
Hamiltonian to produce a time-dependent impurity po- 
tential 



V (t) = V I x + X (t) + ^-,y + Y (t) 



V x (*) 



(6) 



We proceed to switch to the interaction picture through 
the unitary operator Uq — exp (iHot/h) and solve the re- 
maining equation up to second order in time dependent 
perturbation theory obtaining yet another simplified ver- 
sion of von Neumann equation 



ih^ (lAUvWpW ] UlU^ = 0, 



where the time evolution operator is given by 

rt 



(7) 



1 



Vj (si)dsi 

to 
t ps± 

/ V I (s 1 )Vi(s 2 )ds 1 ds 2 . 

to J to 



(8) 



here Vi (t) — UqV (t) U\ is the impurity potential in the 
interaction picture. The formal solution to (0 is given by 
p(t) = W^tfp(t )UU W where p (t ) = p Q = f(H ) 
is the equilibrium density matrix at the initial time to — > 
— oo. 

The density current is proportional to the thermal and 
time average of the velocity operator 



dm [p (t) n] 



(9) 



where S is the surface of the sample, and the limit S — > oo 
is understood. By performing a cyclic permutation in the 
trace we obtain 



J 



-Tr 



(to)««oWnw t « t « t 



(10) 



After lengthy calculations the components of the density 
current is worked out as 

_ ne 2 n Ei - ujcTieijEj 



m 



mm' 

where i,j — x,y and e^j is the antisymmetric tensor 
(e i2 = -e 2 i = 1 and en = e 22 = 0), 



_ N i B\V{q)Y 2 

^mm' — \ u w' \ z > 



and A, 



Smh y qn 

g. t A w / + 2 \ejj\qjUJ^2 
A^+4^ 2 

[Ujq + UJ C (p - (if)] - U} 2 + f] 2 , UJ q 



(12) 



lo x E x + ujyEy, uj x = -nu c (q x + q y TiUj c )/B(l + t 2 uj 2 ), 
LU y = nu> c (-qy + q x TiUj c )/B(l + t 2 uj 2 ) and / M = 
/ (ftw c (// + 1/2)). The matrix elements D 11V are given 
by 




M > M > 



(13) 



where z g = (tfo — iq y )/y/2 and denotes the associ- 

ated Laguerre polynomial. 

Retaining a finite value of the switching parameter r\ 
yields a density of states for the Landau levels with the 
Lorentzian form given in Eq. (|12[) ; it is distorted by the 
electric field through the ui q term. Henceforth we will 
consider rj — Tto c - The differential conductivity tensor is 
calculated from Eq. (fTTj) as <7y = OJi/dEj. Finally the 
differential resistivity tensor is obtained from the inverse 
of the conductivity: that is r»j = cr^ . 

In the limit of small bias and small magnetic field 
the expression for the density current reduces to J x = 
ne 2 TiE x (1 — a) /m where 



a 



2tt 



e ~E F /k B T \V\ 2 N im 
( e E F /k B T + i y 



Shv 2 



(14) 



Hence the quantum scattering time and the inelastic scat- 
tering time can be related by r = rj(l — a) or similarly 
the elastic scattering time is given by r e = Tj(l — a) /a. 
The factor iVi|V| 2 /5T 2 present in the expressions for 
the density current can be estimated from the sample's 
mobility and the inelastic scattering time. 

In a current controled scheme: the longitudinal density 
current is fixed to a constant value Jo while J y should 
vanish. This leads to a set of two implicit equations for 
the density current 



J x (E x ,E y )=J , J y (E x ,Ey) = 0, 



(15) 
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FIG. 3: Electric field as a function of the dc bias J x for 
T = 2if and for fixed magnetic fields ranging from B = 0.5T 
to B — 1.085T. The thin lines indicate that r xx < 0. 
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FIG. 4: Electric field E x as a function of the dc bias J x for 
T = 2K and for magnetic field B = 0.5T. The thin lines 
indicate differential resistivity r xx < 0. The inset shows a 
possible non uniform configuration for the density current. 



where the explicit form of the functions Ji is given in 
Eq. (fTT|) . To obtain the components of the electric field 
E x and E y , we start assigning initial values E x = E Xa 
and E y — E yo that solve these relations in the absence of 
impurities (i.e. using only the first term on the R.H.S. of 
Eq. (|lip ). then the accuracy of the solution is improved 
by a recursive application of Newton's method. 



III. RESULTS 

Fig. [T]shows the differential resistivity r xx — dE x /dJ x 
as a function of the longitudinal dc density current J x for 
a magnetic field B — 0.784T and various values of the 
temperature. We use a sample mobility /i = lOQm^V/s, 
electron density n = 8.2 x 10 15 m -2 and a broadening 
parameter V = 0.04. As the value of the temperature is 
reduced the differential resistance decrease approaching 
zero. We can observe that at low temperature (T < 2K) 
and above a threshold bias current (J x > 0AA/m) the 
differential resistivity becomes negative. Positive values 
for the differential resistance are recovered at higher bias 
or higher temperatures. The strong temperature depen- 
dence observed in this plots, consistent with the experi- 
ments, is originated mainly on the T 2 dependence of the 
inelastic scattering rate. 

The electric field E x is plotted as a function of the lon- 
gitudinal current J x in Fig. O It is important to notice 
that E x differs from the longitudinal voltage by a geo- 
metrical factor. DNRS are observed below T — 4K and 
above the current threshold J x > 0AA/m in the form of 
negative slope curves (see inset of Fig. [2]) in accordance 
with the r xx negative values observed in Fig. [1] Accord- 
ing to Bykov et the stability condition is simply 
expressed as r xx > 0. Thus the regions in Figs. [T]and[2] 
that display a negative differential resistivity are unsta- 
ble, and they should rapidly evolve into ZDRS to insure 
stability. Accordingly in Fig. [1] we should replace the 
NDRS by r xx = and maintain a constant slope in Fig. 
[2] instead of the negative slope. At higher values of J x the 
differential resistivity becomes positive (Fig. [I} as well as 
the longitudinal voltage slope as a result of an increase in 
the impurity scattering prevalent at high electric fields. 
In this regime the large electric field components, neces- 
sary to maintain the strong dc bias and J y = 0, cause 
the impurity terms to strongly participate 2 ^. 

Fig. [3] display a series of plots of E x field as a func- 
tion of the longitudinal density current J x at T = 2K 
for various fixed values of the magnetic field that corre- 
spond to Shubnikov-de Haas oscillations maxima. The 
thin lines indicate negative values of r xx that violate the 
stability condition. As the magnetic field increases the 
width of the electric field plateaus increase and the pos- 
itive slope is recovered for higher onset density currents. 
An isolated plot of the longitudinal electric field E x as a 
function of the dc current J x is shown in Fig. [4] In the 
inset of Fig. [l]we show a nonuniform distribution current 
similar to the one proposed by Bykov et alM. With this 
configuration not only the stability condition r xx > is 
fulfilled but the electric field is uniform throughout the 
sample given that E x = E m i n for 3 X \ and J X 2- The aver- 
age current density J x = [J x \yi + J X 2V2)/{yi + V2) may 
be modulated by varying the sizes y\ and 2/2 of the dif- 
ferent density current domains with the restriction that 
2/1 + J/2 = w. Notice that more complicated schemes with 
more density current modulations also fulfill this condi- 
tions. 
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IV. CONCLUSIONS 

We have presented a model for the nonlinear trans- 
port of a 2DES placed in a strong perpendicular mag- 
netic field. The model is based on the solution of the von 
Neumann equation for 2D damped electrons, subjected 
to arbitrarily strong magnetic and dc electric fields, in 
addition to the weak effects of randomly distributed im- 
purities. This procedures yields a Kubo formula that 
includes the non-linear response with respect to the dc 
electric field. Considering a current controlled scheme, 
we obtain a set of nonlinear self-consistent relations that 
allow us to determine the longitudinal and Hall electric 
fields in terms of the imposed external current. NDRS 
are found in the low temperature (T < 2) and moderate 
bias regime OAA/m < J x < l.QA/m. In low dc bias (low 



electric field regime) the dominant mechanism is the in- 
elastic one. The longitudinal electric field (and voltage) 
recover they positive slope in the high bias (high electric 
field regime) . It is shown that in order to correctly repro- 
duce the main experimental results the inelastic scatter- 
ing rate must obey a T 2 temperature dependence, con- 
sistent with electron-electron Coulomb interaction as the 
dominant inelastic process. 
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